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Abstract: In this paper, we address the atmospheric entry guidance and control 
problem for Mars precision landing. The guidance and control design is based on 
the principle of tracking a reference drag versus velocity profile in the entry flight 
corridor, which is determined by physical constraints of the flight. An integrated 
adaptive /robust control approach to atmospheric entry guidance and control is 
introduced to deal with different uncertainties. 
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1. INTRODUCTION 

To facilitate guidance, navigation, and control de- 
sign for Mars precision landing, the Mars landing 
process is divided into three phases, i.e., atmo- 
spheric entry phase, parachute phase, and terminal 
landing phase. In this study, we only focus on guid- 
ance and control design for the entry phase, i.e., 
the phase from atmospheric entry of Mars orbiter 
to parachute deployment. Feedback aerodynamic 
maneuvering during the entry phase is used to re- 
duce the error to a required level at the parachute 
deployment. Hindering the achievement of this goal 
are uncertainties in the guidance and control sys- 
tems. The main sources of the uncertainties include 
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winds/gusts, air density dispersions, and measure- 
ment and knowledge errors; in addition, some sim- 
plification assumption on mathematical model made 
to facilitate the control design also generates a gap 
the real system and the employed model, which is 
conveniently treated as uncertainty. Basically, the 
above uncertainties are classified as external distur- 
bances and parametric uncertainty. 

In this paper, we will address the guidance and 
control problem from the view point of adaptive 
and robust control design. The following questions 
will be answered to provide an integrated design 
method: given an a priori error bound, does there 
exist any guidance and control law to achieve an 
error level not greater than the given bound under 
all admissible disturbances and vehicle evolution 
parameter changes, and how does one design such 



a control law? In particular, we consider the lon- 
gitudinal landing guidance following McEneaney- 
Mease’s approach [6]. The objective is to achieve 
a certain specified range during entry by control of 
the lift/drag ratio component in position- velocity 
plane, which is achieved by modulation of the Mars 
landing vehicle’s angle of attack and/or bank angle. 
The guidance approach is based on the principle 
of tracking a reference drag-velocity profile in the 
entry flight corridor, which is determined in the 
drag-velocity plane by physical constraints of the 
flight; those constraints include the limitations of 
the thermal protection system surface temperatures, 
normal load factor, dynamic pressure, and equilib- 
rium glide flight (see [2] for details). A guidance 
approach based on drag reference tracking has been 
used in the NASA Space Shuttle re-entry guidance 
law [2,3]. Basing guidance and control law on drag 
is advantageous because drag is essentially directly 
measured by the accelerometer, and the range to be 
flown during entry can be decided by the drag ac- 
celeration profile maintained throughout the entry 
flight and can be predicted using simple analytical 
techniques. 

There exists several guidance and control techniques 
under ideal conditions, i.e., assuming perfect mea- 
surement and no uncertainty [2,6,7], but without 
guarantee for error tolerance level due to distur- 
bances. In this paper, basing on previous guidance 
and control results, we will develop an integrated 
guidance and control law to minimize the impact 
of uncertainties. This 2-part paper is organized as 
follows: in section 2, an uncertain drag dynamical 
model is derived from the basic equations of motion 
with the consideration of atmospheric disturbances 
and some operating parameter variations; the guid- 
ance and control problem for Mars atmospheric 
entry is formulated as a problem of tracking a drag- 
velocity profile. In section 3, the guidance and con- 
trol problem is solved as an adaptive disturbance 
rejection problem where the parametric uncertainty 
and atmospheric disturbances are treated separately 
so as to improve system performance. 


2. UNCERTAIN DYNAMICS OF MARS 
ATMOSPHERIC ENTRY 

The objective of the longitudinal guidance is to 
control the range to be flown during entry, which 
is a function of the drag- velocity profile [2]. The 


main idea for the guidance strategy is to control 
the vehicle to follow a given drag-velocity profile. 

2.1 Equations of Motion 

The equations of translational motion used in this 
paper are developed in a coordinate system with one 
axis oriented along the Mars-relative velocity vector, 
one axis perpendicular to the plane formed by the 
position and Mars relative velocity vectors, and 
the third axis completing the right-hand coordinate 
system. The Mars-relative translational state of the 
landing vehicle is represented by the variables R 
(range), h (altitude), V (velocity relative to the 
Mars), 7 (flight path angle), and xp (heading angle). 
These equations of motion are as follows [10]: 


R = V cosy (1) 

A = U sin 7 (2) 

V = —D - g sin 7 ( 3 ) 

1 V 2 

7 = y (( — 9) cos 7 + Lc os <r) ( 4 ) 


xp = — ( cos 2 7 tan <p sin tp + Lsincr)(5) 

V cos 7 r 

where a is the bank (or roll) angle and <p is vehicle 
Mars-relative latitude, r the distance between the 
lander and the center of Mars, g = p/r 2 with 
H the Martian GM is the Martian gravitational 
acceleration. These equations of motion neglect the 
Coriolis and centripetal accelerations due to the 
Mars’ rotation because they are relatively small. If 
the relative motion of the atmosphere to Mars is 
ignored, the specific drag and lift are given by the 
following: 

D = l p ^ c ° v2; l = \ p ^ Cl y2 ' 

Z m Z m 

where Cd and Cl are drag and lift coefficients, 
respectively; they change with the change of the 
angle of attack and the magnitude of the velocity. 
An exponential atmospheric density p model is used 
in the following. 

2.2 Drag Reference and Range Prediction 

For the Mars atmospheric entry, we employ a drag 
profile in the drag acceleration / Mars-relative ve- 
locity plane. The drag profile is chosen to stay 



within the entry flight corridor due to some physical 
constraints and to minimize the accumulate aerody- 
namic heat load (see [2] for details). A typical drag 
profile, D r = /(V), is shown in Fig. 1, composed 
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Fig. 1. Typical Drag- Velocity Profile 

of a quadratic drag segment, an equilibrium glide 
(linear) segment, and a constant drag segment [2,6]. 
Each of the segments is taken as a general quadratic 
profile: f(V) = Co + c{V + c^V 2 . The objective for 
the entry guidance law is to achieve a desired range 
at a certain specified velocity. In fact, the range can 
be approximated by the following relation: 

V Vf 

R(V) = - f —— - : dV » - f -dV 

J f{V)+g sin 7 J f{V) 

Vo Vo 

at the earlier stage of the flight as the flight path 
angle is small. During the latter part of the entry 
when the landing vehicle is at lower speed, one can 
use an energy variable E = gh + \V 2 to replace 
the velocity variable V in the drag expression [2]. 
Therefore, the range to be flown is determined by 
the drag-velocity (or drag-energy) profile. In the 
following, the control and guidance law is designed 
so that the actual drag will track the given drag- 
velocity profile. 


2.3 Drag Dynamics with Uncertainty 

To facilitate the control design for entry guid- 
ance based on the drag-velocity profile, one should 
directly consider the drag dynamics, instead of 
the dynamical equations considered earlier in (1)- 
(5). We are particularly interested in the impact 


of the atmospheric uncertainties, including surface 
winds/gusts and atmospheric density dispersions. 


Impact of Atmospheric Uncertainties 

The atmospheric uncertainties considered in the 
following include the air density dispersion and the 
winds/gusts. When uncertainties are presented, D 
and L in the equations of motion axe not drag and 
lift, rather they are the aerodynamic forces parallel 
and perpendicular in the plane of symmetry to the 
vehicle velocity vector relative to Mars, respectively. 
The air flow velocity due to winds and gusts can 
be represented as V w = V + AV, where V is the 
velocity of the landing vehicle relative to Mars, 
AV is the disturbance due to the winds. The drag 
and lift coefficients change with the change of the 
angle of attack as well as the magnitude of velocity. 
The influence of the surface winds/gusts, which are 
viewed as disturbances, can also contribute to the 
changes in specific drag and lift coefficients, in which 
case: 

Co = C°D + ACd, C l =C° l +AC l 

where C% and Cj- are the actual drag and lift coef- 
ficients, and A Cd and A Cx are the corresponding 
deviations due to the above treatment. 

The atmospheric density can be modeled as the 
ideal exponential density model (8) plus uncertain 
dispersion A p (see Fig. 2), i.e., it can be well 
represented by the following: 

P = Poe~^+Ap. (7) 
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Fig. 2. Mars Atmospheric Density Dispersion 




The influence of the atmospheric uncertainties con- 
tributes to the changes in specific drag and lift, in 
which case: D — Do 4* A D and L = Lq + A L, where 
Do and Lq are drag and lift without disturbances, 
and AD and A L are the resulting perturbations. 

From the density model (7), one has that 

p — A p _ h 

p - Ap 


D D 2 D D 2 - 

= -^f-^ +2 7 +2 ^- c ’- <5o) - (11) 

From (10), one can solve h; replace the h and 
h defined by ( 9 ) in ( 11 ), one haw the following 
equation: 

= w + gi(D,D,V)u , (12) 


or | = -fa + Sd for some 64 . 

From the drag formulation, one has 

^= 42 ^+^=-|+ 4 +c+M8) 

D p V f Cd H V 

where C := ^ and 8 a = 

c . n AVV-AVV ( A C D C% - A CpC° D 
6d + 2 (V + AV)V + (C° D + AC d )C° d 


Uncertain Drag Dynamical Model 

Now we are ready to derive the drag dynamics 
for control design. The control of the translational 
motion is performed by adjusting u := fjcosa, 
which can be achieved by the modulation of either 
angle of attack or the bank angle. 

In the following, we assume the gravity acceleration 
g = /i/f with r being the reference radius is 
constant. The actual variation of g can be taken 
as a part of the disturbance to be discussed later. 

From the equation (2), (3), and (4), one has 

h = — (V sin 7 ) = V cos 77 + V sin 7 

dt 


for some smooth functions go, gi, and <72 of D, D, 
and V, where 6 is a parameter vector defined by 


0 = 


C 

C 

1 cos 2 7 

r 

1 — cos 7 
sin 2 7 


(13) 


w is the combined disturbances: 

4gHD 2 


w = — 


V(V 2 -4- 2 gH) 


S a 4- D6 a \ 


T/2 I, V* 7* 

= (— - g) - D— - —(1 - -cos 2 7 ) + D cos 7 tl 


With the consideration of atmospheric uncertain- 
ties, from ( 8 ), one has 

(1 + 2 -£j-)h =-»(§+ 27 -c-s.). ( 10 ) 

Take derivatives on both sides, one has 

aH 9HD: A g 2 H sin 2 7 

(1 + 4- 6 y 3 h 4* 4 ^ 


On the other hand, from (3) and (2), one has 

h 

V = -D - 5 sin 7 = -D - g — 

= fo(D, t)y V) + 6 (Dy by V)C + b(Dy D, V)S a (U ) 

= :i{DybyVyCy 8 a )y 

for smooth functions £o> £i > and f of D, D, V', C, and 

<5 0 . 

The above equations give a model for drag dynamics 
with uncertainty. Note that the model (12)-(14) is 
fairly general to cover many other sources of un- 
certainties. In particular, the measurement noises, 
estimation errors, as well as the gravitational ac- 
celeration variations can be conveniently considered 
(9^s part of the combined disturbance w. If we take 
Xi = DyX 2 = Dy and x 3 = V , then complete 
dynamics governing the drag evolution is as follows: 

xi = x 2 

X 2 = -0o(*l>*2 } Z3) “ 0l(Xi,X 2 ,* 3 ,u)0 

+u/ 4 * ^ 2 (xi, x 2 , x 3 )u (15) 

X3 — ^(X1,X 2 ,X 3 ,C, S a ) 

Z — X\ 

In the following, we will use the above dynamical 
equation to design a guidance law. 



3. GUIDANCE AND CONTROL DESIGN 


3.1 Guidance as a Tracking Problem with Disturbance 
Rejection 

As the drag profile is given by the range-velocity 
requirement, to reach certain range requirement, it 
is sufficient for the drag governed by the dynamics 
(12)-(14) to follow the corresponding drag-velocity 
profile in the presence of disturbances. Therefore, 
the guidance problem can be formulated as a track- 
ing problem with disturbance rejection. 

Guidance Problem: Consider the uncertain sys- 
tem (15), suppose 9, € 0, where 0 is a convex closed 
bounded admissible parameter set. Given a drag- 
velocity profile, D r = f{V)\ design a feedback con- 
troller with measurement of D t D, and V, as well as 
the profile information, such that for all admissible 
energy-bounded unknown parameters, D governed 
by (12) satisfies the following requirements: 

• Tracking: if w = 0, 

lim {D(t)-f(V(t))}= 0. 



Fig. 3. Adaptive Disturbance Rejection 

z\ = xi - D r = e, z 2 = x 2 - D r \ 

then equation (15) can be represented as: 

{ z — Az + B\d + B 2 v 
£3 = f (*i + D r , z 2 + br, Z 3 ,C, S a ) (17) 
e = Z! 


• Disturbance Rejection : 

T T 

1 1| D(t) - £> r (t)|| 2 dt<\ 2 j IMOII* dt + e 
0 0 

for given e > 0, and some A > 0. 

In the following section, we will provide an inte- 
grated adaptive/robust control approach to deal 
with the above problem based on the dynamical 
equation (15). The control architecture is illustrated 
in Fig. 3. The techniques of linearization, TL 00 - 
control, and dissipation theory-based adaptive con- 
trol will be used ([5]). 

3.2 Feedback Linearization 

The starting point for this problem is to design a 
controller for the disturbance rejection problem with 
assumption that the parameter 9 is known. Consider 
system (15), take 

Po(Xl, 12,^3) + 5 ?(ll,l 2 ,l 3 )^ + t' /ic\ 

u = t - ; — ^ — ; \ ’ l i0 ' 

(1 + 0 s)g 2 (x u 12,13) 

Z\ 

where v is the new input variable. Let z = , 

where 


where A = 


0 1 
0 0 


,B\ = f?2 = 


0 

1 ’ 


and C = 


[10], The above (error) system is input-output 
linear with linear party state z; we will use a linear 
technique to derive the control v. 


3.3 Hoo -Control Design 

Consider the linear part of the equation (17) with 
state z. Note that the error system is independent 
of the state variable 13 = V and observable from 
output e = z\. Suppose P > 0 is a solution to the 
Riccati equation for the state- feedback TLoo -control 
problem: 

PA + A t P + Pi^BxBf - B 2 B%)P + C T C = 0. 

If A > 1, then P > 0. In this case, the Woo-control 
solution gives feedback controller: 

v = Fz •.= \BjPz. 

with / 0 T (||e(t)l| 2 )dt < A 2 / 0 T (||d(t)|| 2 )<*. One can 
recover the following controller: 

U = K{9-,D r ,Dr,D r ,Xi,Xi,X3) 







Fj_-_Dr +go(xi, Ja,x 3 ) + g?(n, 13,13) 

(1 + 05)52(11, *3.13) 
to achieve the disturbance rejection. 

3.4 Adoptive Control Law 

However, as 9 is not known, the above control law 
can not be directly used. One can use an estimate p 
of 9 instead, and the control law would be 

t i — K (p; Dr » * D r , X \ , , X3 ) ( 18 ) 

where p is given by an update law: 
p = i 7 (x,u). 

Therefore, the guidance and control system has the 
structure as illustrated in Fig. 3. The error equation 
can be equivalently represented as: 

I z = (A + B 2 F)z + B\w 

+ [ , ° J(p-0) 

m 9l(XuX2,Xz } U) ^ 

*3 = Z(zi 4- D r , Z 2 4* £> r ,X3,C,$a) 
e = Zi 

where F is a state feedback matrix defined as F = 
\B%P and P satisfies the following equation: 

P(A + B 2 F) + (A + B 2 F) t P 

+ jpPB\Bi P 4- C T C = 0. 

If = z T Pz, then completion of square implies 

^((^ + B 2 F)z + Bxw) < A a Ml* - ||e|| a (19) 

To design the update law, we take the following 
function as a storage function [5]: 

W(z,p) = V(z) + (p - 9) T Q(p - 9) 

where Q > 0 such that 

max {(p — 9) t Q{ p - 0)} < e. 

p ,*€0 

Therefore, by completion of square, one has 
W(z,p)< A 2 |M| 2 ~ l|el| a + 

2 (jp - 9) T Q (p + Q~ l [0 gf(xi,x 2 ,x 3 ,u) ] Fz) . 


Therefore, if one chooses 

P = *e(P,Q- 1 [0 -g’[(x l ,x 2 ,x 3 ,u) ] Fz), (20) 

where ttq is the vector projection (see [5]), then 

( P ~ 9) T Q {P + Q~ l [ 0 g{ (xi , x 2 , x 3 , «) ] Fz) < 0. 

Thus, we have 

r r 

I ||e|| 2 dt < X 2 j (||tt/|| 2 )dt + t. 

0 0 

Let’s look at the tracking issue if w = 0. We still use 
W{z,p) as the Lyapunov function. Then with the 
above adaptive law, one has W(z,p ) < - ||e|| 2 < 0. 
Therefore, z and p— 0 are bounded. It can be shown 
that / 0 °° ||e|| 2 dt < oo, and e(t)->0 as t-+ oo. 

Therefore, the adaptive control law given by (18)- 
(20) achieves adaptive disturbance rejection. 
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